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Abstract. We review many body calculations of the equation of state of dilute neutron matter in the 
context of effective field theories of the nucleon-nucleon interaction. 
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1 Introduction 

One of the central problems of nuclear physics is to cal- 
culate the properties of nuclear matter starting from the 
two-body scattering data and the binding energies of few 
body bound states PE]. The nuclear matter problem is 
notoriously difficult. Some of the problems that are often 
mentioned are 

— the large short-range repulsive core in the nucleon- 
nucleon interaction 

— the large scattering length in the 1 Sq channel, the small 
binding energy of the deuteron, and the small satura- 
tion density 

— the need to include three (and possibly four) body 
forces 

— the need to include non-nucleonic degrees of freedoms, 
such as isobars, mesons, quarks, etc. 

Ever since the discovery of QCD the classic nuclear 
matter problem has evolved into the broader question of 
how the properties of nuclear matter are related to the 
parameters of the QCD, the QCD scale parameter and 
the masses of the light quarks. 

Over the last couple of year much progress has been 
made in understanding these kinds of questions in the case 
of nuclear two and three-body bound states [3] . Using ef- 
fective field theory methods it was shown that 

— the short range behavior of the nuclear force is not 
observable. Using the renormalization group the short 
distance behavior can be modified without changing 
low energy scattering data and binding energies @E|- 

— effective field theories can accommodate the large scat- 
tering lengths in the nucleon-nucleon system [Sj. The 
scattering lengths depend sensitively on the quark mas- 
ses, see Fig. P and the large value of a( 1 5'o) observed 
in nature appears to be accidental 

— a local three body force is necessary to renormalize the 
two-body force already at leading order. As a conse- 
quence, one cannot predict three-body binding ener- 
gies based on two-body scattering data alone UJ. 



— non-nucleonic degrees of freedom, quark effects, rela- 
tivistic effects etc. can be absorbed in local operators. 

Effective field theories have also achieved remarkable 
quantitative success in describing the available nucleon- 
nucleon scattering data below the pion production thresh- 
old |1UU11| . The long term goal is to achieve a similar 
qualitative and quantitative understanding of the nuclear 
many body problem. 

In this contribution we shall study a simple limiting 
case of the nuclear matter problem. We shall concentrate 
on pure neutron matter at densities significantly below 
nuclear matter saturation density. The neutron-neutron 
scattering length is a nn = — 18 fm and the effective range 
is T nn — 2.8 fm. This means that there is a range of densi- 
ties for which the inter-particle spacing is large compared 
to the effective range but small compared to the scatter- 
ing length. Neutron matter in this regime exhibits inter- 
esting universal properties. We are interested in the limit 
(fci?a nra ) — - > oo and (kpr nn ) — > 0, where fop is the Fermi 
momentum. From dimensional analysis it is clear that the 
energy per particle at zero temperature has to be propor- 
tional to energy per particle of a free Fermi gas at the 
same density 



A ^\AJo q, 5\2mJ 



(1) 



The constant £ is universal, i. e. independent of the de- 
tails of the system. Similar universal constants govern the 
magnitude of the gap in units of the Fermi energy and the 
equation of state at finite temperature. 

Universality also implies that the properties of this sys- 
tem can be studied using atoms rather than nuclei. The 
scattering length of certain fermionic atoms can be tuned 
using Feshbach resonances, see [H| for a review. A small 
negative scattering length corresponds to a weak attrac- 
tive interaction between the atoms. This case is known as 
the BCS limit. As the strength of the interaction increases 
the scattering length becomes larger. It diverges at the 
point where a bound state is formed. The point a = oo is 
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Fig. 1. Quark mass dependence of the scattering length in 
the 1 5o channel. The plot shows a combination of unquenched 
lattice QCD results (red points, from Beane et al. I12| ) and 
chiral extrapolations. The experimental point is shown in pur- 
ple. Two different power counting schemes were employed to 
constrain the quark mass dependence, the BBSvK scheme 6 
113) and the W (Weinberg) scheme 0. 
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Fig. 2. Equation of state of pure neutron matter at T — 4 
MeV from lattice simulations of an effective field theory. Figure 
taken from Lee & Schafer |15|. The curves labeled (fc,fl,f2) 
show results for a free gas on the lattice and in the continuum, 
the curves labeled (bl,b2) show ladder sums, and (sl,s2) are 
numerical results on different lattices. We also compare to the 
variational results of Friedman and Pandharipande (FP). 



called the unitarity limit, since the scattering cross section 
saturates the s-wave unitarity bound a — 4ir/k 2 . On the 
other side of the resonance the scattering length is posi- 
tive. In the BEC limit the interaction is strongly attractive 
and the fermions form deeply bound molecules. 



2 Numerical Calculations 

The calculation of the dimensionless quantity £ is a non- 
perturbative problem. In this section we shall describe 
an approach based on lattice field theory methods. The 
physics of the unitarity limit is captured by an effective 
lagrangian of point-like fermions interacting via a short- 
range interaction. The lagrangian is 



C = ipU id, 



The standard strategy for dealing with the four-fermion 
interaction is to use a Hubbard-Stratonovich transforma- 
tion. The partition function can be written as |15) 



Z = J DsDcDc* exp [-S] 



(3) 



where s is the Hubbard-Stratonovich field and c is a Grass- 
mann field. S is a discretized euclidean action 



'Cj(n)ci(n, + 6) 



_ e v'-Coa tS (r l )+-%^( 1 _ 6 /j) c |( n ) c .( n ) 

h ^2 [ c i( n ) c t( n + h) + c*(n)ci(n - l s ) 



(4) 



Here i labels spin and n labels lattice sites. Spatial and 
temporal unit vectors are denoted by l s and 0, respectively. 
The temporal and spatial lattice spacings are b T and b. The 
dimensionless chemical potential is given by ft = fib T . We 
define a t as the ratio of the temporal and spatial lattice 
spacings and h — at /(2m). Note that for Co < the 
action is real and standard Monte Carlo simulations are 
possible. 

The four-fermion coupling is fixed by computing the 
sum of all two-particle bubbles where on the lattice. Sche- 
matically, 

mil v ^ 1 

4^ ~~ Co + 2 t—~EJ 
v H 



(5) 



where the sum runs over discrete momenta on the lattice 
and E p is the lattice dispersion relation. A detailed dis- 
cussion of the lattice regularized scattering amplitude can 
be found in [T?)lll7lll5| . For a given scattering length a the 
four-fermion coupling is a function of the lattice spacing. 
The continuum limit correspond to taking the temporal 
and spatial lattice spacings b T , b to zero 



b T [i — > 0, 



bn 1/3 -> 0, 



(G) 



where /i is the chemical potential, n is the density and 
an 1 / 3 is fixed. We performed numerical simulations at 
non-zero temperature and concluded that £ = (0.09 — 
0.42). Lee studied canonical T = simulations and ob- 
tained £ = 0.25 |18|. Green Function Monte Carlo calcu- 
lations give £ = 0.44 an( i finite temperature lattice 
simulations have been extrapolated to T — to yield sim- 
ilar results [HH1II - 

Lattice results for the equation of state of dilute neu- 
tron matter at T = 4 MeV are shown in Fig. [21 For com- 
parison, we show variational results obtained by Friedman 
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and Pandharipande using a phenomenological potential 
|22| . We observe that the lattice calculations agree very 
well with the variational result. The pressure is very sim- 
ilar that of non-interacting neutrons scaled by a factor 
~ 1/2. The lattice calculation can be extended to higher 
densities by including explicit pionic degrees of freedom 
in the effective lagrangian [231 - ^ n this case a mild sign 
problem returns, but at T / this sign problem can be 
handled with standard methods. 





Fig. 3. The particie iadder diagrams 
the leading order contribution to the 
the large d limit. 



4 Large d expansion 



shown in this figure give 
ground state energy in 



3 Analytical Approaches: Large N expansion 

It is clearly desirable to find a systematic analytical ap- 
proach to the dilute Fermi liquid in the unitarity limit. 
Various possibilities have been considered, such as an ex- 
pansion in the number of fermion species |24U25| or the 
number of spatial dimensions 26 27,28 22j . 

We begin with a brief description of the large N ap- 
proach. The physics of the large N limit depends on the 
symmetries of the interaction. One possibility is a SU (N) 
symmetric interaction |24| 



(7) 



where / = 1,...,N is a flavor label. A smooth large 
N limit is achieved by keeping C$N — c constant as 
N — ► oo. The large N limit is most easily studied by in- 
troducing a Hubbard Stratonovich held a coupled to the 
density iptipf. The leading contribution to the free energy 
comes from the free fermion term and the mean field con- 
tribution, both of which scale as N. Subleading 1/N cor- 
rections arise from particle-hole ring diagrams. The prob- 
lem is that at any hxed order in the large N expansion 
the free energy diverges as the scattering length is taken 
to infinity. 

This problem is related to the fact that particle lad- 
ders need to be summed in the unitarity limit. This can 
be achieved by by considering a Sp(2N) symmetric inter- 
action of the form 



C 



Go 



(8) 



where J = (02) <8 • • • ® (02) and Go A = g constant 
as A — > 00. This interaction can be bosonized using a 
difermion field <P = (ip fj j ' 9 tp g ) /A . At large A the leading 
contribution corresponds to the mean field BCS approx- 
imation. The thermodynamic potential in the unitarity 
limit is 



Q 



-A 



d 3 p 
(2tt) 3 



e; + <P 2 



m<P 2 



(9) 



with e p = E p — /1. This function can be minimized nu- 
merically. We find $0 = 1.16// and £ = 0.59. Nikolic and 
Sachdev studied 1/N corrections near T c 25 . These ef- 
fects are not small. They find, for example, (/j,/T)(T c ) = 
1.50 + 2.79/A + 0(1/N 2 ). 



Steele suggested that the many body problem of non- 
relativistic fermions near the unitarity limit can be studied 
using an expansion in 1/d, where d is the number of spatial 
dimensions |26j . The main idea is that phase space factors 
associated with hole lines are suppressed as d — > 00 so 
that the leading order contribution comes from 2-particle 
ladders, and higher order corrections correspond to the 
hole line expansion of Bethe and Brueckner. 

Consider the effective lagrangian in equ. (J2J). We first 
study pcrturbative corrections to the ground state energy 
in d spatial dimensions. The leading order correction to 
the energy per particle is 



A 



n d c k d f 2 M 

(2w) d 



( 



V 2M 



(10) 



This expression indicates that the large d limit should be 
taken in such a way that 



A 



Q d C k d F - 2 M 
d(2n) d 



const. 



(11) 



In the following we wish to study whether this limit is 
smooth even if the theory is non-perturbative. Consider 
the in medium two-particle scattering amplitude in d spa- 
tial dimensions. The result is 



d d q 

{2n) d ~W 



fvac(k) 



2(2vr) 



fpp(K, s). 



(12) 

The theta function = 9(ki — kp)9(k2 — kp) with fci.2 = 
P/2 ± k requires both fermion momenta to be above the 
Fermi surface. The first term on the RHS is the vacuum 
contribution. In dimensional regularization the vacuum 
term is purely imaginary and does not contribute to the 
ground state energy. The second term is the medium con- 
tribution which depends on the scaled relative momentum 
k = k/kp and center-of-mass momentum s = P/(2kp). 
In the large d limit we find 



/pp( S ! K ) — ^fpp( S l K ) 



1 



•(0), 



1 



(13) 



which implies that all two-particle ladder diagrams are 
of the same order, see Fig. The sum of all ladder di- 
agrams can be calculated by noting that, except for the 
logarithmic (BCS) singularity at s = 0, k — 1, the particle- 
particle bubble is a smooth function of the kinematic vari- 
ables s and k. Hole-hole phase space, on the other hand, 
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is strongly peaked at s 



- k — 1/V2 in the large d limit. 
4/d • (1 + 0{l/d))). The ladder 
sum is a simple geometric series and |27) 



We find that f p d) P (s,R) 



-r = U + 



1 - 2A 



"'3 



\,2M / 



(14) 



where A is the coupling constant defined in cqu. f| 1 1 p . We 
observe that if the strong coupling limit A — > 00 is taken 
after the limit d — > 00 the universal parameter £ is given 
by 1/2. We have also studied the role of pairing in the 
large d limit. The pairing gap is 



„ 2e^E F 
A = : exp 



1 

dX 



1 + 



(15) 



There is no exponential suppression in d — > 00 limit, but 
AjEp is down by a power of 1/d. As a consequence the 
pairing energy is sub-leading compared to the result in 
equ. (TH|) . 



5 Epsilon expansion near four dimensions 

Nussinov & Nussinov observed that the fermion many 
body system in the unitarity limit reduces to a free Fermi 
gas near d = 2 spatial dimensions, and to a free Bose gas 
near d = 4 [2Hj . Their argument was based on the behavior 
of the two-body wave function as the binding energy goes 
to zero. For d = 2 it is well known that the limit of zero 
binding energy corresponds to an arbitrarily weak poten- 
tial. In d = 4 the two-body wave function at a = 00 has a 
1 /r 2 behavior and the normalization is concentrated near 
the origin. This suggests the many body system is equiv- 
alent to a gas of non-interacting bosons. 

A systematic expansion based on the observation of 
Nussinov & Nussinov was studied by Nishida and Son |291 

. In this section we shall explain their approach. We be- 
gin by restating the argument of Nussinov & Nussinov in 
the effective field theory language. In dimensional regu- 
larization a — > 00 corresponds to Co — > 00. The fermion- 
fermion scattering amplitude is given by 



\ m J 



r ( 1 -|) {~po + E p /2~i5Y- 



(16) 



where 8 — > 0+. As a function of d the Gamma function has 
poles at d = 2, 4, . . . and the scattering amplitude vanishes 
at these points. Near d — 2 the scattering amplitude is 
energy and momentum independent. For d = 4-ewe find 



A{p ,p) 



8ir 2 e 



m 2 po — Ep/2 + id 



+ 0(e 2 ). (17) 



We observe that at leading order in e the scattering am- 
plitude looks like the propagator of a boson with mass 
2m. The boson-fermion coupling is g 2 = (8tt 2 e)/m 2 and 
vanishes as e — > 0. This suggests that we can set up a per- 
turbative expansion involving fermions of mass m weakly 
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Fig. 4. Leading order contributions to the ground state energy 
in the e = 4 — d expansion. Solid lines are fermion propaga- 
tors, dashed lines are boson propagators, and the cross is an 
insertion of the chemical potential. 



coupled to bosons of mass 2m. In the unitarity limit the 
Hubbard-Stratonovich transformed lagrangian reads 



L = & 



id + <r 3 — 



2m 



(18) 

where 'P = (ip^,tp^) T is a two-component Nambu-Gorkov 
field, o"j are Pauli matrices acting in the Nambu-Gorkov 
space and a± — (a± ± ia%)j2. In the superfluid phase <j> 
acquires an expectation value. We write 

a A , VSn 2 e ( m4> V /4 

4> = 0o + g<p, g = -7T- > ( 19 

m \ 2tt J 

where (fig = (cf>). The scale M 2 = m4>o/(2ir) was intro- 
duced in order to have a correctly normalized boson field. 
The scale parameter is arbitrary, but this particular choice 
simplifies some of the loop integrals. In order to get a well 
defined perturbative expansion we add and subtract a ki- 
netic term for the boson field to the lagrangian. We include 
the kinetic term in the free part of the lagrangian 



C = ^ 



id + <7 3 - h 4> Q ((T + + <7_) 



2m 
V 2 

4m 



(p. 



(20) 



The interacting part is 

Li = g (W^a+Wip 



h.c) 



9 ' ld ° + 4^ ' ^~ 



(21) 



Note that the interacting part generates self energy correc- 
tions to the boson propagator which, by virtue of equ. (|17|) . 
cancel against the kinetic term of boson field. We have also 
included the chemical potential term in Li. This is moti- 
vated by the fact that near d — 4 the system reduces to a 
non-interacting Bose gas and /1 — > 0. We will count \i as a 
quantity of 0(e). 

The Feynman rules are quite simple. The fermion and 
boson propagators follow from equ. I|20(l and the fermion- 
boson vertices are iga . Insertions of the chemical poten- 
tial are ifia 3 . Both g 2 and \i are corrections of order e. In 
order to verify that the e expansion is well defined we have 
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to check that higher order diagrams do not generate pow- 
ers of 1/e. Studying the superficial degree of divergence of 
diagrams one can show that there are only a finite number 
of one-loop diagrams that generate 1/e terms. 

The leading order diagrams that contribute to the ef- 
fective potential are shown in Fig. 0] The first diagram is 
the free fermion loop which is 0(1). The second diagram 
is the (i insertion which is O(l) because the loop diagram 
is divergent in d = 4. The sum of these two diagrams is 



(27T) 



E l 




(22) 



The integral can be computed analytically. Expanding to 
first order in e = 4 — d we get 



Vn = 



1 



7-3( 7 + log(2)) 
6 



1 l-2( 7 -log(2)) c 



m<p 
2ir 



d/2 



(23) 



Nishida and Son also computed the two-loop contribution 
shown in Fig. 0] The result is 



V 2 



-Ce 



A \ d / 2 



2- 



(24) 



where C ~ 0.14424. We can now determine the minimum 
of the effective potential. We find 



2// 



[l + (3C-l + log(2))e + 0(e 2 )] 



(25) 



The value of V — V\ + V2 at (f> determines the pressure 
and n — dPj dp gives the density. We find 



1 

n = — 

e 



l-i(2 7 -l-log(2))e 



m(f)o 
2tt 



d/2 



(26) 



We can compare this result with the density of a free Fermi 
gas in d dimensions. This equation determines the relation 
between ep = fcf,/(2m) and the density. We get 



2tt 




U 










-)] 


m 




v2" 





2/<l 



(27) 



We determine 6f for the interacting gas by inserting n 
from cqu. I|26l) into cqu. I|27|l . The universal parameter is 
£ = fi/eF ■ We find 

£ = ie 3/2 + 4:e 5/2 Iog(e) - 0.025e 5/2 + . . . = 0.475, (28) 
2 16 



O OO 
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Fig. 5. Leading order contributions to the ground state energy 
in the e = d — 2 expansion. 



6 Epsilon expansion near two dimensions 

Near two spatial dimensions the scattering amplitude in 
the unitarity limit vanishes linearly in e — d — 2 



2tt 

A(p ,p)=i—e + O(e 2 ). 
m 



(29) 



The coefficient of e is momentum end energy independent. 
This means that we can set up a perturbative expansion 
with an effective four- fermion coupling g 2 = 27re/m. This 
expansion is very similar to the perturbative (kpa) expan- 
sion studied by Huang, Lee and Yang |32II33| . see Fig. [5] 
but it is not restricted to the weak coupling limit. To 0(e) 
the effective potential is given by 

V0 + V1 = -P free -^p 2 + 0(e 2 ), (30) 

where P/ ree is the pressure of free fermions expanded to 
0(e) and the density is given by 



(31) 



From the total pressure we can compute the density and 
Fermi energy as in the previous section. The universal 
parameter £ is given by 



£= 1-1+ 0{f 2 ) =0 (e = l). 



(32) 



Similar to the perturbative expansion pairing is exponen- 
tially suppressed in the e expansion. The pairing gap is 

m 

to = ^exp(-iVl + 0(e)), (33) 

which corresponds to the perturbative result of Gorkov 
and Melik-Barkhudarov [3^ . Equation (|32(l shows that the 
e expansion is poorly convergent. However, the e expan- 
sion is useful in improving the convergence of the e = d — 4 
expansion, and in connecting the perturbative (kpa) ex- 
pansion with the physics of the unitarity limit. 



where we have set e = 1. The calculation has been ex- 
tended to 0(e 7/2 ) by Arnold et al. pH]- Unfortunately, 
the next term is very large and it appears necessary to 
combine the expansion in 4 — e dimensions with a 2 + e 
expansion in order to extract useful results. 



7 Outlook 

In this contribution we focused on an idealized systems 
of neutrons at very low density. The obvious question is 
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to what extent these methods can be extended to nuclear 
systems near saturation density. 

The lattice simulations can easily be extended to in- 
clude finite range effects, explicit pions, and isospin. Some 
of these refinements will cause a sign problem in the simu- 
lation, but in most cases the sign problem can be handled 
with standard methods. A significant amount of work will 
be required in order to reduce discretization errors to the 
point where the interactions are quantitatively reliable all 
the way up to momenta on the order of the Fermi momen- 
tum in nuclear matter. 

The large N, large d, or epsilon expansions are easily 
extended to interactions with a finite scattering length and 
a finite effective range |25p35lf3"n1 . There is also no obvious 
obstacle to including explicit pion degrees of freedom. It 
will be interesting to extend these methods to systems of 
protons and neutrons. In this case three-body forces have 
to be included. The central question is whether saturation 
can be achieved, and whether effective field theories pro- 
vide a qualitative understanding of the Coester line \X7\ . 

We should also note that the many body physics that 
governs the equation of state of nuclear matter near sat- 
uration density may well be simpler than the physics of 
the unitarity limit. Effective range corrections suppress 
the two-body scattering amplitude and nuclear matter is 
more perturbative than dilute neutron matter. As a con- 
sequence, perturbative calculations using soft potentials 
or effective interactions adjusted to nuclear matter prop- 
erties may well be reliable [HEIEHI- 
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